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The imbalance of dijets produced in hadronic collisions has been used to extract the average 
transverse momentum of partons inside the hadrons. In this paper we discuss new contributions 
to the dijet imbalance that could complicate or even hamper this extraction. They are due to 
polarization of initial state partons inside unpolarized hadrons that can arise in the presence of 
nonzero parton transverse momentum. Transversely polarized quarks and linearly polarized gluons 
produce specific azimuthal dependences of the two jets that in principle are not suppressed. Their 
effects cannot be isolated just by looking at the angular deviation from the back-to-back situation, 
rather they enter jet broadening observables. In this way they directly affect the extraction of the 
average transverse momentum of unpolarized partons that is thought to be extracted. We discuss 
appropriately weighted cross sections to isolate the additional contributions. 

PACS numbers: 12.38.-t; 13.85.Ni; 13.88.+C 



I. INTRODUCTION 



Event shape observables have been widely studied for various reasons. In e + e~-annihilation, observables such as 
the thrust and jet broadening have been studied primarily to extract a s (M z ), cf. for instance Refs. [l], H, H, 0, H, Q 
for theoretical studies and Refs. @, H, H, EH EH EH EH EH Ell for experimental studies. In the center-of-mass system 
(cms) of the e + e~ collisions at lowest order in ev s , the produced quark-antiquark pair is exactly back-to-back leading 
for two-jet events to a thrust T equal to unity. Gluon radiation, i.e. order a s corrections, gives rise to nonzero 1 — T 
and also to nonzero jet broadening. In the perturbative regime these observables can be used to extract a s , which has 
been done recently at next-to-next- to- leading order [lj| E3, Ell ■ The results compare very well with those obtained 
by other means of extraction. In the nonpcrturbative regime, hadronization will also lead to nonzero event or jet 
shapes. This is characterized by a mean transverse momentum leading in general to a contribution suppressed 

by a power of the large scale, the cms energy Q. For example, in the nonperturbative regime 1 — T oc (k±)/Q. It 
has been suggested that this contains universal information on a s in the infrared regime. We refer to Ref. 19] for a 
review on this topic. 

Event and jet shapes have also been studied in hadronic collisions. Compared to e + e~-annihilation here the 
additional complication of initial parton transverse momenta arises. Another difference is that instead of the thrust 
axis, it is common to use the transverse thrust axis nt, which is the axis in the transverse plane having maximum 
transverse energy flow. The corresponding transverse thrust is defined as [2(1 : 



= 1 



where p l T is the transverse momentum of the outgoing hadron i, Et — \Pt \ is the total transverse energy (neglecting 
masses) and the transverse thrust axis is the transverse unit vector n t that maximizes T t . Here we use the notation 
of Ref. {20j, where also the jet broadening variable Qt is defined as: 

n 

Qt = Y,\PT xn t\- (2) 
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An experimental investigation of the average Q t as a function of E? in pp collisions has been reported in Ref. [2111 . 
Higher order perturbative corrections to the transverse thrust and jet broadening are discussed in e.g. Refs. [22ll23j|. 

Assuming collinear factorization and ignoring broadening from hadronization, Q t will be zero for 2 — > 2 partonic 
subprocesses and only sensitive to 2 — > 3 processes, like for e + e~-annihilation except that there are more subprocesses 
to consider in hadronic collisions. Extraction of a s in hadronic collisions [3, [2f| is however complicated due to the 
presence of parton transverse momenta and the transverse momentum distribution of hadrons inside the jet. The 
former effect one can minimize by considering events with at least three pronounced jets, which means considering only 
large values of Qt, whereas the latter effect could be minimized by considering Q t for jets, instead of hadrons. In fact, 
the quantity Q t for two-jet events, where i now denotes the i-th jet and n = 2, has been used to study and extract the 
average parton transverse momentum. This has been done for instance in Refs. [26l. l27l l28l. l29l. l30l. l3ll. |32[ l33l. |3~S, [35| . 
As can be seen from those results, the average parton transverse momentum extracted from the data increases with 
energy (^/s) and is in general much too large to be attributable to "intrinsic" transverse momentum. This is a 
consequence of soft parton radiation, similar to what happens for the transverse momentum distribution of vector 
boson production in hadronic processes [3|| [37], [H, He| (see also the instructive discussion in Ref. [13]). Resummation 
of soft radiation effectively broadens the transverse momentum dependence of the parton distributions, increasingly 
so with increasing center of mass energy. 

In this paper we point out that besides initial parton transverse momentum and soft parton radiation, there are 
additional contributions to Q t , even for the simplest two-jet case. These are contributions due to the transverse 
polarization of quarks and the linear polarization of gluons inside the initial unpolarized hadrons. These contributions 
can arise for nonzero initial parton transverse momentum. We will show how these effects contribute to Qt and discuss 
that besides complicating the extraction of the average parton transverse momentum from Q t , they may even hamper 
that extraction altogether depending on their magnitude. 

In Ref. [2(| collinear factorization was assumed, making the observable {Qt) only sensitive to 2 — * 3 subprocesses. In 
reality collinear factorization is not always applicable, due to the partonic transverse momentum effects. In a simple 
picture of a Gaussian distribution of intrinsic parton momentum p t , the average value (pt) can be extracted from {Qt) , 
but in fact, no factorization theorem has been established for two-jet or two-hadron production in pp or pp collisions for 
observables that are sensitive to parton transverse momenta. To make matters worse, in the framework of transverse 
momentum dependent parton distribution functions, nowadays commonly referred to as TMDs, it even seems that 
factorization cannot be established for this particular type of process when taking into account nontrivial effects of 
gauge links [4l|, [42|, HH, [H, [45[. This would cast doubt on any conclusion drawn from (Q t ) in hadronic collisions, 
except for large Qt where collinear factorization can be applied. But even if factorization will work out in some as 
yet unknown way, the additional contributions from spin dependent TMDs may complicate matters considerably. 
Schematically this can be seen as follows. 

Consider the process h\h,2 — > jijzX, where ji stands for produced jet i. In the plane transverse to the collision 
axis, 6(j> denotes the deviation of the (azimuthal) angle between the two jets from tt, i.e. 8<j) = <pj 1 — <f>j 2 — ~k. It is 
sometimes referred to as the dijet imbalance. Let us consider only 2 — > 2 subprocesses. In collinear factorization the 
6<p dependence of the cross section will then only receive a contribution at 8cj> = 0. Allowing for parton transverse 
momentum in the initial hadrons leads to a smearing of the 8<f> distribution. For the idealized case of equal jet 
transverse momenta (both equal to Et/2) the differential cross section takes the form: 

A(Ql)+B(Ql)Ql + C{Q 2 M, (3) 



dE T d8<t> 



where Q t = Et\ s'm(8(f>/2)\ is equal to the absolute value of the transverse momentum of the two-jet system. A, B, C 
are functions of Q\ , which do not need to vanish at Q\ = 0. The terms B and C appear from spin effects inside the 
initial hadrons hi, for which expressions will be presented in this paper. In general these spin-dependent contributions 
are not suppressed by powers of 1/ Et, also not when arising from polarized gluons as claimed in Ref. [46| . A result 
for B has recently been obtained in [46[ following a calculation similar to the one for pp — > jjX presented in [47l |. 
This contribution arises from the quark TMD hf q (48|, which represents the distribution of transversely polarized 
quarks inside an unpolarized hadron. The new result in this paper is the contribution from h^ 9 [4^ |. the distribution 
of linearly polarized gluons inside an unpolarized hadron, which gives rise to C. Upon ignoring these spin effects, only 
the term A remains and the average Qt value in that case will indeed be directly related to the average transverse 
momentum that is thought to be extracted in Refs. [2(| H3, [H, H^, [3(1 HH, [H, [H, [H, [35J. Our results in principle 
cast doubt on whether the actual value of (kf) has been extracted in those cases. In practice, it all depends on the 
magnitude of B and C. We will present a simple Gaussian model to illustrate the generic shape of the modification 
of the dijet imbalance distribution by B and C terms. 

The paper is organized as follows. First we will present the calculation and expressions for the cross section in Eq. 
([3]) , assuming factorization in terms of transverse momentum dependent correlators and ignoring the possible effects 
from gauge links. We will actually discuss the more general case in which the two jet transverse momenta are not 
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equal, but differ by a small amount w.r.t. Ef. In that case the angular dependence is more involved than given in Eq. 
© , even upon expansion in the small transverse momentum difference of the two jets with respect to their sum. We 
will first express the cross section in terms of the individual jet momenta through their sum and difference (section HT1 
in particular Eq. (|16[) ) and subsequently in terms of the sum and difference of the lengths of the jet momenta in order 
to arrive at the dijet imbalance distribution expressed in more standard variables (section IIII| in particular Eq. I|51|) ). 
In section HVl we discuss angular-projected asymmetries, such as (cos Sip) and the ones that can be used to extract 
B and C. After that we consider the consequences of nonzero functions for the jet broadening quantity Qt, in 
particular for the averages (Qt) and (Qt). Finally (section IVlj) we briefly address the open issues of factorization 
(breaking) and color flow dependence upon inclusion of gauge links. We end with conclusions and two appendices, 
one on relations among various variables in the transverse plane and one on photon-jet production that completes the 
treatment given in Ref . [47j . 



II. THEORETICAL FRAMEWORK: CALCULATION OF THE CROSS SECTION 



We consider the process 

h 1 (P 1 )+h 2 (P 2 )^)et(K 1 )+iet{K 2 )+X, (4) 

where the four-momenta of the particles are given within brackets, and the jet-jet pair in the final state is almost 
back-to-back in the plane perpendicular to the direction of the incoming hadrons. Along the lines of Ref. [47|], we will 
instead of collinear factorization consider a generalized factorization scheme taking into account partonic transverse 
momenta. We make a lightcone decomposition of the two incoming hadronic momenta in terms of the light-like 
Sudakov vectors n + and n_, satisfying = ri 2 . = and n+-n- = 1: 

M 2 M 2 
Pf = P+< + ^+nt , and P% = —^n^ + P^nt . (5) 



2Pf 2P 
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The partonic momenta (pi, p 2 ) can be expressed in terms of the lightcone momentum fractions x 2 ) and the 
intrinsic transverse momenta (pit, P2t), as follows 



rf = x 1 P 1 + < + ^-§21^+rf T , and p^ P -fJ^L^ ++ x 2 P 2 n^ + p^ T . (6) 



In general n+ and n_ will define the lightcone components of every vector a as a* = a ■ n^, while perpendicular 
vectors a± will always refer to the components of a orthogonal to both incoming hadronic momenta, Pi and P 2 . 
Therefore in Eq. ([6|), if we neglect hadron masses, p^ T = and p^ T = Pax- We denote with s the total energy 
squared in the hadronic cms frame, s = (Pi + P 2 ) 2 = E 2 ms , and with iji the pseudo-rapidities of the outgoing partons, 
i.e. 7]i — — In (tan(i6>;)), Oi being the polar angles of the outgoing partons in the same frame. Finally, we introduce 
the partonic Mandelstam variables 

s = {pi+P2?, t=(pi-Ki) 2 , u={pi-K 2 f, (7) 

which satisfy the relations 

t 1 i & 

-- = v = t i > and -- = l-y. (8 

Following Refs. and we assume that at sufficiently high energies the hadronic cross section factorizes in a 
soft parton correlator for each observed hadron and a hard part: 

$ a (xi,p 1T ) ® <f> b {x 2 ,P 2 T) ® \H ab ^ cd (pi,p 2 ,Ki,K 2 )\ 2 . (9) 

This form assumes the simplest possible factorization omitting any gauge link dependence in the correlators, which 
can modify or even break the factorization (see section IVII for a discussion of these open issues) . 

In Eq. ([9]) the sum runs over all the incoming and outgoing partons taking part in the reaction. The convolutions 
(g> indicate the appropriate traces over Dirac indices and \H\ 2 is the hard partonic squared amplitude. The parton 
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correlators are defined on the lightfront LF (£-n = 0, with n = n_ for parton 1 and n = n + for parton 2); they 
describe the hadron — > parton transitions and can be parameterized in terms of transverse momentum dependent 
(TMD) distribution functions. In particular, the quark content of an unpolarized hadron is at leading twist (omitting 
gauge links) described by the correlator [48] 



(2tt) 3 



f< (p\mm\p)\ 



(10) 



where f1(x,p^) is the unpolarized quark distribution, which integrated over p T gives the familiar lightcone momen- 
tum distribution f1(x). The time-reversal (T) odd function hi <1 (x,p 2 n ) is interpreted as the quark transverse spin 
distribution in an unpolarized hadron [48| . Analogously, for an antiquark, 



1 



<pi m m \p) j LF = - m x ,p 2 T ) r + 9 >,pD 



2M 



(11) 



The gluon correlator (omitting gauge links) is given by [4 

•d(£.P) d 2 £ T 



(p-n) 2 
1 

2~t: 



(2tt) 3 



e *P< (p| Tr[F^(0)F^(O] \P)\ 



LF 



ft fl(x,p%) + ) h^i^pl) \ , 



V M 2 



2M< 



with gj, being a transverse tensor defined as 



(12) 



(13) 



The function /f(x,p^) represents the unpolarized gluon distribution, while the T-even function hf 9 (x,p 2 ^) is the 
distribution of linearly polarized gluons in an unpolarized hadron. 

In order to derive an expression for the cross section in terms of parton distributions, we insert the parametrizations 
in Eqs. (fTU|) , (fTTj) and (fl"2"l) of the TMD correlators into Eq. Furthermore, utilizing the decompositions of the 
parton momenta in Eq. ©, the <5-function in Eq. © can be rewritten as 



8\pi+P2-K 1 -K 2 ) = -6[ |Jfi 



61 x 2 -^ F (\K 1± \e-^ +\K 2± \e^) 



x5 2 (p 1T +p 2T -K 1± -K 2± ) 



(14) 



with corrections of order 0(1/ s). After integration over x\ and x 2 , from the first two (^-functions on the r.h.s. of Eq. 
(fT4")h one obtains 



\K 1± \e^+\K 2 



o'/2 



X 2 - 



\K 1± \e-^+\K 2 



->ii 



(15) 



which relates the partonic momentum fractions x\, x 2 to the rapidities and the transverse momenta of the jets. These 
basic tree- level relations will be used in our treatment. We will not consider several other effects that need to be 
accounted for in practice such as the actually used jet definition and higher order corrections that affect the above 
relations and cause additional smearing. 

The hadronic cross section can be written in the form 



( ^ cr /ii?i2^jetjetX 



„2 r 



drndr] 2 d 2 K 1± d 2 K 2± S K 2 , 



MfA) + b (Qt)Qt cos2(0 t - (/)_!_) + C(q^)q% cos4(0 T - <j> ± ) 



(16) 



where q T = K±± + K 2 j_ and K± = (K±±_ — K 2 ±)/2. The sum momentum q T is useful as an in principle accessible 
experimental observable momentum which in our calculations via the delta function in Eq. (|14|) is related to intrinsic 
transverse momenta, q T = p 1T +p 2 T- We denote with (f>T and <p±_ the azimuthal angles of q T and K±, respectively. 
Besides q^, the terms A, B and C depend on other kinematic variables often not explicitly indicated, namely y, x\, 
x 2l and contain convolutions of the various parton distributions. These are discussed separately in the following three 
subsections, where explicit expressions for them can be found, calculated at leading order (LO) in perturbative QCD. In 
deriving these expressions we will often employ the approximation |q> T | <C ~ |-^2±| ~ I-^jJ which is applicable 

in the situation in which the two jets are almost back-to-back in the transverse plane. However, in deriving Eq. (|16[) 
we must be particularly careful with the angular dependence, because approximations in the angular dependence 
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that boil down to approximating <j>i sa <\> 2 + it (such as in Eq. (21) of Ref. [47]) will of course not give the proper 
dependence of the dijet imbalance angle 8<f> = 4>\ — (j) 2 + n. In Eq. (|16p the combination <fij< — 4>± appears, which will 
allow to isolate the terms B and C by q^-weighted integration over q T (cf. section ITVT) . However, in order to arrive 
at the 5<p distribution discussed in the introduction, it is more convenient to express the cross section in terms of 
the combination 4>t — where 4>j is the average jet direction angle, i.e. <pj = (<Pi + 4>2 — t)/2 with <pi and <j) 2 the 
azimuthal angles of the two outgoing jets in the transverse plane. In the present case where \K±\ ^> (|J£"ix| — |-K"2x|), 
it holds that <f)± w <f>j allowing the two angles to be identified to good approximation for all values of 6<j> (cf. Eq. 
(|A8j) ). In the limiting case when |-Kij_| = |-K"2x|, the angles <f)± and (f>j exactly coincide and the T and _L directions 
are orthogonal, so we have exactly cos2(</>t — </>x) = — 1 (note that this will lead to Eq. with a minus sign in 
front of B, but that is of course only a matter of definition) and cos4(0t — <f>±) = 1- This implies that all angular 
dependence then resides in q\, which is in that case solely depends on the off-collinearity of the jets through the dijet 
imbalance angle 8<fi (discussed in section UlTj) . 



A. Angular independent part of the cross section 



The term A in Eq. (|16p is the angular independent part of the cross section and is given by the sum of several 
contributions A ah ^ cd coming from the partonic subprocesses ab — > cd underlying the reaction h\ h 2 — > jet jet X: 



A(y, Xl ,x 2 ,q 2 T ) = A ab ^ cd {y,x 1 ,x 2 ,q 2 T ). 



(17) 



with a,..., d = q, q', q, q' , g. We denote with q and q' two quarks having different flavors, and similar notation holds 
for the antiquarks. Furthermore, the following convolutions of unpolarized parton distributions are defined 



T a \x 1 ,X 2 ,q 2 T ) 



d p 1T d p 2T S (p 1T + p 2T - q T )fi{x 1 ,p 1T )f 1 {x 2 ,p 2T ) , 



(18) 



where a sum over all (anti)quark flavors is understood. Our results for the terms A ah ~* cd in Eq. (TIT]) are listed below, 
starting from the ones corresponding to the (anti)quark induced processes, 



A qq 



= a(y)F qq (x 1 ,x 2 ,q 2 T )+a(l-y)F qq (x 1 ,x 2 ,q 2 T ). 



(19) 



j^qq -^qq 
J^qq-^qq 
j^qq—>qq 
J^qq^q'q' 



a{y)T qq (x u x 2 , q%) + o(l - y)T q g (x u x 2 , q%) . 



N 2 -l 
2N 2 



y(i - y) 



\ + {l-y) 2 1 + y 2 



1 



(1-y) 2 Ny(l-y) 



f qq ( Xl ,x 2 ,q 2 T ) 



b(y)T qq ( Xl ,x 2 ,q 2 T ) + 6(1 - y)J* q (xi,X2,<£) , 



N 2 -l 
2N 2 



2/(1 -y) [y 2 + {i-y? 



F qq ( Xl ,x 2 ,q 2 T ) + T qq (xi,x 2 , q%) 



A qq - 



■<J'J 



N z -U 2 2 1 \ y 2 + (1 - y) 2 

y + (i - y) - 



N \" ' v a ' N 2 
with N being the number of colors and 

N 2 -I 1 + (1 - yf 

a (y^ = 2N 2 ^ ~ v > y ' b (y) = a (y) + 

Analogously, from the gluon induced processes, one has: 

A*™ = c(y)^( Xu x 2 ,q 2 T )+c(l-y)^ q (x 1 ,x 2: q 2 T ), 



a— = iNNii:^,,,^, 



F qq ~(x U X 2 , q 2 T ) + F qq (xi,x 2 ,q 2 T ) 



N 2 -l 
2N 2 



2/(1 - y) 



2 ^n ,2 2(1-^ 
y +(i-2/) 



N 2 - 1 2/(1 - y) 



(20) 
(21) 
(22) 
(23) 

(24) 

(25) 

(26) 
(27) 
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A"- 



■qq _ 



where 



N ( 2 . , 2 1 \ y 2 + (1 - y) 2 , 2x 



(28) 



1 + (1 - yf 



1 + (1 - y) 2 y 

y iV 2 



(29) 



Agreement is found between the results given in the present subsection and the explicit expressions of the partonic 
cross sections published, for example, in [5l], \E^, [H, [SH, HH . However, with respect to Ref. ,461 we find agreement 
with the expression for the unpolarized qq production subprocess, but not for the qq production subprocesses. In 
particular, we find differences as compared with their Eqs. (24) and (33). 



B. The cos2(4>t — 4>±) angular distribution of the dijet 



In Ref. [46] it is shown that the subprocesses qq —* qq and qq — > qq contribute not only to the angular independent 
part of the cross section, according to Eqs. (|2Tj) and ([22]) . but also to an azimuthal asymmetry of the dijet arising 
from the product of two T-odd functions, h^ q hf q or h^ q h^ q . Such an asymmetry is similar to the one calculated in 
the Drell-Yan [56| and in the photon-jet production |47[ processes. We refer to [47j for the details of the derivation 
and present here only our final results. In analogy to Eq. (|17p . we write 



B{y 1 x 1 ,x 2 ,q 2 T ) = B ° 6_>cd (f> a; l) a; 2,9r)) 



with 



B" 



N 2 -l 
N 3 



y(l-y)W 1 (x 1 ,x 2 ,q 2 T ), 



= d(y)n^(x u x 2 ,q 2 T )+d(l-y)H^(x u x 2 ,q 2 T ), 



and 



d(y) 



N 2 



N 3 



■ y(l - y) 2 (1 + Ny) 



(30) 

(31) 
(32) 

(33) 



The following convolution of (transversely polarized) quark and antiquark distributions has been introduced 
q 2 T U qq {x ll x 2 ,q 2 T ) = M ^ M2 J d 2 p 1T d 2 p 2T 5 2 (p 1T +p 2T -q T ) 

(2{h ■ p 1T )(h ■ p 2T ) - (p lT -p 2T )j h^ q {x 1 ,p 2 1T )hf q (x 2 ,pl T ) , 



flavors 1 



(34) 



with h = 9t/|9t|) an d a similar definition holds for 7i qq upon replacement of q — > q in Eq. (f34|) . The small-q T 
behavior of Tt is regular provided the integrations over pj, h 1 9 (x, p 2 -,) con verg e. In addition to Eqs. (|3"Tj) and (p?2"j) . 
we find that the subprocesses qq — ► gg and qq — > q'q', not considered in [46j, also show a cos2(0t — 4>±_) angular 
dependence, leading respectively to 



99^99 — 



and 



B 



2399^9 9 



N 2 -l 
N 



N 2 



(±-yY-jp)y(±-y) 



n qS (x 1 ,X2,q 2 T )+n qg (x 1 ,x 2 ,q 2 T ) 



N 2 



n qq {x 1 ,x 2 ,q 2 T ) + n qq {x 1 ,x 2l q 2 T ) 



(35) 



(36) 



Agreement is found between the results given in the present subsection and the explicit expressions of the polarized 
partonic cross sections published in [54], HH, H3, HH • For the polarized qq production subprocess we find agreement 
with Ref. 46], but again not for the qq production subprocesses (in particular, we find a difference compared to their 
Eq. (26))/ 
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C. The cos4(</> T 



angular distribution of the dijet 



The cos4(0t — 4>±) angular distribution of the dijet is related to the presence of linearly polarized gluons in 
unpolarizcd hadrons. This being a new result of the present paper, its derivation will be discussed in some more 
detail. The gluon-gluon induced part of the reaction under study, to lowest order in pQCD, is described in terms of 
the partonic two-to-two subprocesses 



g{Pi) + aijPi) -> g(Ki) + g(K 2 ) , and 
The corresponding cross sections are given by 



g{pi) + 9{P2) -» q{K x ) + q(K 2 ) 



(37) 



da 99 ^ 99 



dr]i dr\ 2 d 2 K ±d 2 q T S K A 



A"^"{y, xi, X2,gx) + / d 2 p 1T d 2 p 2T S 2 (p 1T +p 2T - q T ) 



N 2 



N 2 — I 



y(l - y)(l - y(l - y)) T"(p 1Tl p 2T , K 1U K 21 _) 



(38) 



and 



da 99 ^ 99 



drji drj 2 d 2 K± d 2 q T S K 2 



A"^ 99 (y,x 1 ,x 2 ,q? r )- / d 2 p 1T d 2 p 2T 6 2 {p 1T +p 2T - q T ) 



N 2/(1 -y) 



N 2 - 1 



y 2 + (l-yf-l^)p^ {Prr .p 2l .K l± .K 2 _) 



(39) 



where 



T"( PlT ,p 2T ,K 1± ,K 2± ) = 



- PitP2T + -77T (K-lx ■ K 2 ±)(p 1T ■ p 2T ) - {Kxx ■ Pit)( k 2± ■ P 2T ) 



K' 



(Ki± -p 2T )(K 2L -p 1T ) 



M 2 M 2 h i 9 Oi » Pit) h i 9 ( x 2 , P 2 2T ) 



(40) 



The functions A"~* 99 and A"~^ 99 , given in Eqs. (|27|) and (|28j) . contain the convolution of unpolarized gluon distri- 
bution functions T" defined in Eq. (|18p. In order to show that the two cross sections in Eqs. ([55)1 and (|3"5)) can be 
written in the same form as Eq. (|16[) . we introduce the functions 



q\,I 3g (xi,x 2 ,q%) = M l M 2 J d 2 p 1T d 2 p 2T 6 2 (p 1T +p 2T - q T ) 

x (2(h ■ p 1T )(h ■ p 2T ) - (p 1T ■ p 2T )^j hl 9 {xi,p\ T )h^ 9 {x 2 ,p 2 2T ) , 

where again h = 5t/|<7t|, and 

q^C"(x 1 ,x 2 ,q 2 ^) = M 2 M 2 J d 2 p 1T d 2 p 2T 5 2 (p 1T + p 2T - q T )p 2 1T pl T hf 9 { Xl ,p\ T )h^ 9 {x 2 ,p 2 2T ) . 

The small-q T behavior of 1 and £ are regular provided the integrations over pj, h^ 9 {x,p\) converge. Hence we have 
J d 2 p lT d 2 p 2T S 2 {p 1T + p 2T - q T )V" 

n) r 

2 ( S il§im _ § ij s lm + pmgjl^gg _ jgg^ 



(41) 



(42) 



4 E 



K {l K 3} K {1 K m} 



{i j} \ / {I m} 



q% 



{21 99 - C 99 ) 



= q A T cos2(20 T - 0i - 4> 2 )(2J 99 -C 99 ). 



(43) 



The difference between the angular dependence cos2(2(/)x — 4>i — §2) — cos4(0r — 4>j) and cos4(0x — 4>±) is °f order 
q 2 T jK\ (cf. Appendix©. Substituting Eq. 05} into Eqs. J38J) and ([39]), and defining da 99 = da 99 ^ 99 + da 99 ^ 99 , we 
finally obtain 

2 

A 99 (y,x 1 ,x 2 ,q 2 T )+C 99 (y,x 1 ,x 2 ,q 2 T )q^ cos4(0 T - 4> ± ) 



da 99 



dm dr)2 d 2 K ± d 2 q T S K 2 



(44) 



where A" = A 99 ^ 99 + A 99 ^^, C" = C"^ 9g + C 99 ^™, with 
N 2 



and 



C'~' 



Qgg^qq 



N 2 - 1 



y(l-y)(l-y(l-y)) 



2I 99 (x 1 ,x 2 ,q' T )-jC 99 (x 1 ,x 2 ,q^) 



N y(l-y) i .> 

'v 2 + (i-y) 2 - 



N 2 - 1 



1 

N 2 



21"(xi,X2,q 2 T )-C 99 (xi,x 2 ,q 2 T ) 



(45) 



(46) 



It turns out that the two subprocesses gg — > gg and gg — > qq are the only ones that determine the cos4(0y — fa\_) 
dependence of the cross section. Therefore in Eq. (fTB]) 



which, together with Eqs. l(iT ]) -((i2 ll and Eqs. 1(35 1 -(|i6 ll . leads to 



C 



N 



N 2 - 1 



y(i - v) 



i 



N(l-y(l-y))- 1 [y 2 + (l- y y - — 



2I"(x u x 2 ,q T ) - £"(x 1 ,X2,q 2 T ) 



(47) 



(48) 



showing how the azimuthal asymmetry under investigation is related to the T-even, spin and transverse momentum 
dependent parton distribution function hf 9 (x,p T ). 



III. DUET IMBALANCE DISTRIBUTIONS 



In this section we study the cross section for the process h\ h,2 — > jet jet X in terms of the total transverse energy 
Et and the dijet imbalance 5(f> = fa — fa — it, which are the kinematic variables commonly used in the experiments. 
The dijet imbalance angle describes the deviation of the two jets from a back-to-back configuration (see Fig. [5] in 
Appendix |A"|) . 

The transverse energy is the sum of the transverse energies of the two jets, Et = |-Ki_l| + I-K'axl) an d the difference 
is defined as AA'i = |-K"i±| — |-K~2±|- In our basic expression for the cross section in Eq. ([To]) we have traded Ki± 
and K 2 ± for q T and K±, but we can also trade the variables (|1Cll|, |-K"2_l|) for (Et, AK±) and (fa, fa) for (fa, 
6 fa). We find in the back-to-back approximation 



q% 



4Ki 



AK 2 cos 2 | ^ 



EUos<l 6 4 



Ei sin 2 -f 



AK 2 sin 2 | ^ 



AK 2 



Ei 



Ei sin 2 -f 



(49) 
(50) 



In the first expression we cannot drop the term proportional to AK 2 ^ because it is not a good approximation for 
5<j) m 0, which is most relevant. Note also that this implies q\ > AK\, i.e. AK\ sets a lower bound on the q\ values 
probed, which may be very relevant if the functions A, B, C are steeply falling functions with increasing q T . 
The cross section in Eq. (fT7j)) rewritten yields 



^./ii/^^jetjetX 



A(q 2 T ) + B(q 2 T ) q 2 T cos2(<?(> T 



dr]idrj2dETdAK±dfad5<j) 2s 
with q T given in the unapproximated first part of Eq. (|49p and 

q\ cos2(4>t - fa) = AK\ cos 2 



■C(q T )q T cosA(fa - fa) 



q T cos4(0t — fa) 



El sin 2 



E T sin 2 



— - AK 2 cos 2 



1 2 



E T AK\ sin 2 (8 fa) 



(51) 

(52) 
(53) 



In this way we have arrived at an expression that is amenable to phenomenological studies, approximating \Kn 
\K2±] ~ Et/2 only in places where the difference is negligible for all values of 5(f>. 
For AK± = we obtain: 



c £ a .h 1 h 2 ^jetjet: 



dr]idrj2dETdAK±d(f>jd5<j) 2s 



A(q 2 T )~B(q 2 T )q 2 T + C(q 2 T )q T 



(54) 
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with in that case exactly q\ = Ej, sin 2 (Sip/2), and in essence recovering Eq. §5§ (the sign in front of B is just a matter 
of definition). 

To illustrate the effect of nonzero B and C terms, we will make a Gaussian Ansatz for these functions of q\. We 
will take: 

c>2 n4 c>6 

A(q 2 T ) = ^exp(-q 2 T R 2 A ), B(4) = ^exp(-4i?|) ) C(q 2 T ) = -f- exp(-q 2 T R 2 c ), (55) 
normalized such that 

J d 2 q T A{q 2 T ) = 1, j d 2 q T q 2 T B(q 2 T ) = 1/c, J d 2 q T q 4 T C(q 2 T ) = 1/c 2 . (56) 

Fig. Q] shows a plot of the cross section in Eq. (f5Tj) as a function of 6<f> for the arbitrary, but perhaps realistic choices 
l-Kijl = 30 GeV, \K 2 ±\ = 31 GeV, R A = 0.5 GeV~ 2 , R B = 2R A ,R C = 3R A and c = 3. For smaller AK± the 




FIG. 1: An illustration of the effect of sizeable B and C terms on the 8<f> distribution of the cross section in Eq. (|51[) . 

shoulders become more pronounced, but already for AK± = 2 GeV the shoulders are hardly distinguishable anymore. 
In general, B and C have to be significant in size and broad enough to generate an observable effect, i.e. for the 5<fi 
distribution to deviate visibly from a Gaussian distribution. 

Although the B and C terms were not considered before in experimental analyses of dijet imbalance measurements 
in hadronic collisions, there is experimental data available that has some bearing on the size of B compared to A. It 
comes from the measurement of the violation of the Lam- Tung relation in the Drell-Yan process. As shown in Ref. 
[561 ] . this violation n is given by the ratio q 2 T W q /T qq (Eq. §M§ divided by the angular averaged result in Eq. (|18|) ) 
but with the sums over flavors weighted with a factor e 2 , the quark charge squared. This has the effect of emphasizing 
the contribution from up quarks. In the present two-jet production case, the ratio q\ B/A in the midrapidity region 
(?7i ~ ?7 2 ~ 0) and for large N can be approximated by q\ H qq jT qq . The size of n in Drell-Yan may thus be expected 
to give some indication of the size of q\ B/A. The violation of the Lam- Tung relation in Drell-Yan has recently been 
measured in pp and pd collisions [59] ]. It is consistent with no violation, but with sizeable errors. Small violation 
would be in line with the expectation that for antiquarks inside a proton is considerably smaller than for quarks. 
For pp one however expects a large violation, as observed in np collisions [13, Q . So the effect of a nonzero h, 
for quarks may be mostly relevant for jet broadening studies in pp [2ll ] . 



IV. WEIGHTED CROSS SECTIONS 



Apart from the fact that nonzero functions for quarks and gluons modify the 8<f> distribution and hence affect 
the extraction of the average initial parton transverse momentum from this dijet imbalance distribution, it would in 
principle be of interest to extract these functions themselves from it. Therefore, the question arises whether one can 
project out the B and C terms separately. In Ref. [46| this is discussed for B only, but there are some problems 
with the proposed method. It was suggested that (cos 8(f>), i.e. the cross section integrated over 5<j) weighted with an 
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additional factor of cos 6<fi, projects out a contribution from hf q exclusively 1 . However, our result in Eq. (|5ip shows 
that (cos 5(j)) does not project out B nor a part of B exclusively, not even in the idealized case when |-K"ix| = l-^2_i_|> 
as can be seen from Eq. (|54p . To see the appropriate weighting, we return to the form in Eq. I|16p and note that B is 
projected out by 



^ ( jh 1 h 2 — >jet jet X y q/2 

(cos2(0 T -^)). oos2^ T -M ^^^^^ -^TBiy^Wr) ^ (57) 

2 



Integrating over the length of q T gives with possible inclusion of additional weighting with powers of 



„2 \ M r [n 1 \ J_/ii/i 2 -»jctjct X 



in which we get for M = 1 the factorized result 

flavors 



(59) 



in the B contributions. For M = 0, the expression does not deconvolute. In that case usable, but model dependent 
expressions may be obtained by making a Gaussian Ansatz for the transverse momentum shape of h\. For the type 
of convolution that appears in B this has been done in the literature, see for instance Ref. (63j . 

Next, we will analyze in some more detail the weighted asymmetry that projects out C . A measurement of the 
weighted cross section 

/dd> ^0-^1^2— *jet jet X ^ ^2 

■W cos4( ^ - ^ dm dm ¥k1 d 2 q T = 2 4 ^'^ 4) ' (60) 

with C given in Eq. (f48|) . would give access to the linearly polarized gluon distribution of a hadron. After integration 
over the length of q T with possible inclusion of additional weighting with q\, we obtain 

r ( a 2 \ M f ( a 2 \ M dcr /iift 2 ^jctjctjsf 

(l&j <-4(^-^)>=/A T (^j -^-^)^-^-^-. (61) 

In this case we get for M — 2 the deconvoluted result 

tt | dq 2 T (j^j 4 - O") = 96C (2) (^i)C (2) (^), (62) 

in the C contributions. In order to evaluate the integral in Eq. (l60|) without weights or study the explicitly in- 
dependence, one can employ a Gaussian model for h^ 9 , of which the easiest choice has a factorized x and p T 
dependence, that is, neglecting the dependence on the factorization scale, 

hi 9 (x,p 2 T ) = ^ hi 9 (x) e- R * p T , (63) 

- / (#r)" = T^yr ^(«), (64) 

where i?^ is a size parameter related to the average partonic p\ by the relation R\ = l/(p§,). For incoming 
(anti)protons, i? p = Rp = R, so one has 

J d 2 q T q i T (21" ~ C") = j^f2 ^ d 2 q T d 2 p 1T d 2 p 2T S 2 (p 1T + p 2T - q T ) 2^{h ■ p lT ){h ■ p 2T ) - {p lT ■ p 2T )^j 



PitP\t 



R*_ 

TT 2 



ht 9 { Xl )hl; 9 {x 2 )e- R2 (PiT+P2T) . ( 6 5) 



1 In Ref. [4ull actually (P^/M 2 cos Stfi) was considered, where Px ~ |^i±| ~ I -Rax I j despite the fact that iC2X was integrated over. The 
factor P?/M 2 artificially enhances the weighted asymmetry if not divided by (P?/Af 2 ). 



11 



Using the p 2 T integration to eliminate the delta function in Eq. (|65p and shifting the integration variable p 1T 
Pit = Pit ~ \Qt^ one arrives at 



d 2 q T q% (22" - C 9 



1 



I? 4 



M 2 M 2 16 7T 
1 



d 2 q T dp'i T q T e 



-R 2 (2pg 



54 h i 9 ( x i) h i 9 0^) • 



Substituting Eq. 
result 



2M2M| i?4 "i ^' • (66) 

into Eq. (foT))) shows that for a Gaussian shape one finds for the unweighted average the final 



7T y dq| <cos4(0t - X )) = ^j- 1 y(l 



y) 



1 



iV(l-y(l- 2/ ))-- r (^ + (l-^-_ 



1 



(67) 



V. JET BROADENING 



In our almost back-to-back jet situation, the jet-direction j coincides with the transverse thrust axis and the jet 
broadening variable Qt defined in Eq. |2]) is given by 



Qt = Et 



\q T \\sm(<pr 



|q T ||sin(0 T - <£_!_) I , 



(68) 



for which we refer to Eq. (|A6|) in Appendix E] and one needs to use Eqs. (|A7p and (IA8|) to check the validity of 
the approximation. Using this expression we can now turn to the evaluation of the average jet broadening (Qt) as a 
function of \K 1 I, 



(Qt) a 



2- 



d 2 q T Qt {\q T \,<h,<t>±) 



dg-hihi— >jet jet 



A" 



(69) 



The differential cross section in the integrand is obtained from Eq. (| 16[) and contains, besides the well-known, angular 
independent term A, also the terms B (due to the transverse polarization of quarks and antiquarks in the colliding 
hadrons) and C (related to the linear polarization of gluons). The following integrals, 



d4>T I sm(<pT 



2tv 



,2tt j 

/ d(j) T I sm((j) T - (j)±)\ cos4(0 T - (j> ± ) = - — 
Jo 15 



ot I sm((pT 



|cos2( 



-'-5 



are all different from zero, meaning that the A, B and C terms contribute to (Q t 



(Qt) cx J d 2 q T \q T \ 



A(q 2 T )- l -B(q 2 T )q 2 T -^C(q 2 T )q T 



In order to calculate (Q 2 ), one needs to evaluate the integrals 

r 2ir 



ot sin 



t>T sin 



&y sm 



0±)cos2(0 T - </>J_) 



7T 

"2 



cos4(0 t - ^_l) = . 



which show that only the terms A and B enter in the estimate of (Q 2 ): 



(Qt) OC /" d 2 q T q 2 T 



MQt) ~ ^B{q 2 T )q 2 T 



(70) 
(71) 
(72) 

(73) 

(74) 
(75) 
(76) 

(77) 
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VI. COLOR FLOW DEPENDENCE AND FACTORIZATION 

In our treatment in this paper we have simply convoluted the quark and gluon correlators with the hard partonic 
cross sections, without worrying about possible nontrivial effects arising from the gauge link structure in these cor- 
relators. The proper gauge invariant definitions of TMDs as well as collinear correlators involve nonlocal operators 
containing path-ordered exponentials, the gauge links. The gauge link is the result of resumming all gluons with po- 
larizations along the momentum of a particular hadron into the soft parts. In the case of TMDs the path of the gauge 
links generally depends on the process. The path dependence disappears after integration over transverse momenta. 
In the collinear correlators, one can usually choose a gauge that makes the gauge link unity, but the same procedure 
for TMDs can leave transverse pieces that are situated at lightcone infinity. These links can have physical effects, for 
instance in single transverse spin asymmetries that arise from the Sivers effect, which is described by a T-odd TMD. 
The Sivers asymmetries in semi-inclusive deep inelastic scattering and the Drell-Yan process are predicted to differ 
by a sign as a consequence of the gauge links [64| . 

In the more complicated processes h\h,2 — > 7/ jet + jet + X the single spin asymmetries involving the Sivers func- 
tion f&BL [6^ | come from correlators with more complex paths in the gauge links. This causes deviations that are 
more involved than a simple sign change with respect to e.g. semi-inclusive deep inelastic scattering. But also in this 
case, calculable process-dependent "color flow" factors can be obtained which may be different for each hard partonic 
subprocess. In this way they allow for the calculation of particular weighted cross sections in dijet production, re- 
sulting in a small asymmetry [6(1 (6^, [68| , as also shown by the data [69[ . However, claims of possible factorization 
breaking have been put forward for this process [iH |44| and this remains an open question. 

For observables involving a product of two T-odd TMDs, such as the one discussed in the present paper, the 
situation is less clear. For cos2</> asymmetries in Drell-Yan [56[ and h\h,2 — ► TjetX the effects of nontrivial gauge 
links were included in Ref. [4tJ following the methods outlined in Refs. [U [42|, HH, [7(|. In both cases the color 
flow factor obtained was +1. However, since the methods used were developed for observables involving a single 
non-contracted transverse momentum px for a T-odd TMD for one of the hadrons in the process, the extension to 
cases in which non-contracted transverse momenta of partons in two different hadrons are involved certainly needs 
careful study. The j»T-dependence for in the correlators for gluons, moreover, has a rank two tensor structure in 
the non-contracted transverse momentum, although it is T-even. For the present case of dijet production (for which 
nontrivial color flow factors were presented in Ref. [46[), which is necessarily more complicated and for which doubts 
about factorization have been put forward, at this stage we do not include any color flow factors. Since we have 
presented the expressions for each partonic subprocess separately, it is possible to include the correct factors at a later 
stage, once they have been firmly established. If factorization cannot be proven for the process of interest, however, 
this not only implies that the functions cannot be extracted but neither can in that case (p T ) be obtained. 

VII. SUMMARY AND CONCLUSIONS 

In this paper we study the effects of transverse momenta of the initial state hadrons in hadronic dijet production. 
The transverse momentum produces an imbalance in the dijets in the transverse plane. In the usual treatments 
the effects are attributed to gluon radiation and to the transverse momentum dependence of the unpolarized quark 
distributions. We look at the effects of two additional TMD functions that enter in the scattering of unpolarized 
hadrons, the distribution of transversely polarized quarks {h^ q ) and the distribution of linearly polarized gluons 
(h^ 9 ). They produce specific azimuthal dependences of the two jets, that are not suppressed a priori. The effects 
cannot be isolated by only looking at the angular deviation from the back-to-back situation, but depend on the jet 
transverse energy and the contributions to it of the two jets. We have discussed appropriate weighting to isolate the 
specific additional contributions. We also pointed out their effect on the jet broadening quantities (Qt) and (Qf), 
which we considered for the simplest two-jet case, but the conclusion that functions contribute to them also affects 
the more general cases containing a sum over hadrons. This in principle complicates the extraction of the average 
initial parton transverse momentum from the jet broadening, but possibly even hampers it altogether if factorization 
of the (sufficiently sizeable) spin dependent contributions indeed turns out to be broken. 
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APPENDIX A: TRANSVERSE PLANE VARIABLES 



In the transverse plane we have the two jet momenta K%± and K 2 i_ defining azimuthal angles 0i and 02. From 
them one can construct the sum and difference angles, 



<j>i = (0i + 2 - tt)/2, 

(50 — 01 — 02 — IT. 

The sum and difference of the transverse energies of the two jets, |.Ki_l| and |-K"ij_|, define 

E T = \K 1± \ + \K 2 ±\ 
AK x = \K tl ]-\K 21 ] 

One can use d 2 K\±_ d 2 K 2 ± = h {Ej. — AKl) dEx dAK± d(f>j d5(f> for the phase space or go to the sum and dif- 

q 



(Al) 
(A2) 

(A3) 
(A4) 




FIG. 2: The transverse plane is defined as orthogonal with respect to the two incoming hadrons. The jet direction (j) is defined 
as (j>j — (0i + 4>2 — 7r)/2. The momenta q T = K 1 ± + K 2 ± and K± = (Ki± ~ K 2 ±)/2 define the azimuthal angles 0t and <j>±- 

ference momenta and their angles as shown in Fig. [21 In that case one has d 2 Ki± d 2 K 2 ± = d 2 K± d 2 q T = 
\q T \ \K±\d\q T \ d\K±\ dcpT d(f>±. We have the following exact relations 



\q T \ cos(0t — 4>j) = AK± cos((50/2), 






(A5) 


\q T \ sin(0 T -(f)j) = E T sin(<50/2), 






(A6) 


2 \K ± \ cos(0_l - 4> 3 ) = E T cos((50/2), 






(A7) 


2 \K X \ sin(0 ± - fa) = AK± sin(<50/2), 






(A8) 


q% = AKl coa 2 (S(t)/2) + E\ sin 2 ((50/2) = AK\ + (E% - 


AKl) 


sin 2 ((50/2), 


(A9) 


4 K\ = E% cos 2 ((50/2) + AK\ sin 2 (60/2) =E\- {E% - 


AKl) 


sin 2 ((50/2), 


(A10) 


2|q T | \K ± \ cos(0 T - ± ) = E T AK±, 






(All) 


2\q T \ \K±\ sin(0 T - <f> ± ) = (E% - AK\)/2 sin(<50). 






(A12) 



We note that the order of the momenta is \q T \ ~ AK± ~ M (hadronic scale) while Et ~ 2 \K±\ ~ y/s, so we see 
from Eq. (| A6|) that (50 ~ M/y/s and from Eq. (|A8j) that 0j_ — <j>j ~ M 2 /s. Further useful relations are 

(A13) 
(A14) 
(A15) 



|q T | 2 cos2(0 T - 0±) w |q T | 2 cos2(0 T - 4>j) = AKl cos 2 ((50/2) - £§. sin 2 ((50/2), 



|q T | 2 sin2(0 T - 



\q T \ 2 sin2(0 T - 0^-) = E T AK± sm(<50), 



|q T | 4 cos4(0 T - ± ) ps |q T | 4 cos4(0 T - 3 ) ps AiT*_ + sin 4 ((50/2) - 6 £f sin 2 ((50/2). 



APPENDIX B: PHOTON-JET PRODUCTION 

In this appendix we include expressions for the terms A and _B for the photon-jet production case, because in Ref. 
[47| we only considered approximate angular dependence. Similarly to Eq. (|16|) . one can write 



cva s 



drj 7 drjjd 2 K 1 ±d 2 Kj± sK\ 



A(q 2 T ) + B(q 2 T )q 2 T cos 2(0 T 



(Bl) 



with 



A(y, a*, a*, <#) = *4 99 ^ 9 + A^ 9 , Bfo, n, x 2 , q 2 T ) = B«~ 



■~ y 



(B2) 
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By comparison with Eqs. (15), (16), and (19) in Ref. [471 ] . we find the following expressions 



h(y)^ a (x 1 ,x 2 ,q 2 T ) + h(l - y)^(x 1 ,x 2 ,q 2 T ) 



with 



h(y) = -(l-y)(l+y 2 ), 



N 2 -l 



and 



N 2 - 1 



f™( Xl ,x 2 ,q 2 T ) +T™( Xl ,x 2 ,q 2 T ) 



Bq ^ g = 2 __ y{1 _ y) J2e 2 q 



H q Hx 1 ,x 2 ,q 2 T ) + H^(x 1 ,x 2 ,q 2 T ) 



(B3) 

(B4) 
(B5) 

(B6) 



in agreement with the results in Refs. [54], 157J, [56 
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